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^'■ 
^ : ABSTRACT 

I We demonstrate that nucleosynthesis in rapid, high-entropy expansions of 

^ ' proton-rich matter from high temperature and density can result in a wider va- 

Q"^ ' riety of abundance patterns than heretofore appreciated. In particular, such 

rO I expansions can produce iron-group nuclides, p-process nuclei, or even heavy, 

^ ■ neutron-rich isotopes. Such diversity arises because the nucleosynthesis enters 

^ ' a little explored regime in which the free nucleons are not in equilibrium with 

O ■ the abundant ^He. This allows nuclei significantly heavier than iron to form in 

1-^ ■ the presence of abundant free nucleons early in the expansion. As the temper- 

I ' ature drops, nucleons increasingly assemble into ^He and heavier nuclei. If the 

;_i ■ assembly is efficient, the resulting depletion of free neutrons allows disintegration 

^ ■ flows to drive nuclei back down to iron and nickel. If this assembly is inefficient, 

IJ • then the large abundance of free nucleons prevents the disintegration flows and 

V~j ■ leaves a distribution of heavy nuclei after reaction freezeout. For cases in be- 

j_j ■ tween, an intermediate abundance distribution, enriched in p-process isotopes, is 



c^ 



frozen out. These last expansions may contribute to the solar system's supply of 
the p-process nuclides if mildly proton-rich, high-entropy matter is ejected from 
proto-neutron stars winds or other astrophysical sites. Also signiflcant is the fact 
that, because the nucleosynthesis is primary, the signature of this nucleosynthesis 
may be evident in metal poor stars. 



1. Introduction 

Apart from numerous important studies of explosive hydrogen burning in matter ac- 
creted onto white dwarfs or neutron stars and early attempts to understand the p-process 
nuclei in terms of hydrogen burning in supernovae (Burbidge et al. 1957; Audouze & Truran 
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1975), astrophysicists have focused relatively little attention on explosive nucleosynthesis 
in proton-rich stellar environments. There are two principal reasons for this. First, after 
hydrogen burning, mainline stellar evolution proceeds at conditions of equal numbers of neu- 
trons and protons, or, because of weak interactions, at a slight degree of neutron richness. 
Accordingly, if the star subsequently explodes, the resulting nucleosynthesis typically oc- 
curs in a neutron-rich environment. Only if the exploding matter has not finished hydrogen 
burning will the nucleosynthesis typically be proton-rich and proceed in a series of proton 
capture and /S^-decay reactions know as the rp-process (Wallace & Woosley 1981). Such 
nucleosynthesis occurs when matter accreted onto the surfaces of white dwarfs or neutron 
stars explodes under degenerate conditions, giving rise to novae or X-ray bursts. 

The second reason is that the nucleosynthesis is thought to be already relatively well 
understood. The details of the rp-process have been well-studied (e.g, Schatz et al. 2001) 
since it was first delineated. For proton-rich matter that achieves higher temperatures, it is 
generally imagined that the proton-rich nucleosynthesis would freezeout from an equilibrium 
in which iron-group isotopes would dominate the abundances. The underlying equilibrium 
abundance distribution might be slightly modified by proton captures at the end of the 
burning. 

The purpose of this letter is to show that explosive nucleosynthesis in proton-rich envi- 
ronments can be much more complex than previously thought. While it is indeed the case 
that the final abundances for many conditions are dominated by iron-group isotopes, for suf- 
ficiently fast, high-entropy expansions of proton-rich matter, isotopes considerably heavier 
than iron can form. These other distributions of nuclei can include interesting quantities 
of light and heavy p-process isotopes or even neutron-rich species usually ascribed to the 
r-process. 

We have not identified an astrophysical site in which such fast, high-entropy, proton- 
rich expansions may occur. The setting envisioned, however, is a neutrino-driven wind 
from a high-mass proto-neutron star early in its epoch of Kelvin-Helmholtz cooling by neu- 
trino emission. Some calculations find wind entropies s/ks of order 100-200 and expansion 
timescales (which we here define as the radial expansion timescale, r = r/f , where r is the 
radial coordinate and v is the radial velocity of a parcel of wind matter) as short as a few 
milliseconds (e.g, Thompson et al. 2001). The Ye-, that is, the electron-to-baryon ratio, in 
these winds is set by the interaction of neutrinos and antineutrinos with free nucleons. While 
the antineutrino spectrum from the neutron star is considerably harder than the neutrino 
spectrum at late times, which tends to make the wind neutron rich, it is possible that the 
two spectra are more nearly equal earlier. This could allow for proton-rich (that is, Y^ > 0.5) 
matter since the lower mass of the proton would favor the reaction Ue + n ^ p + e~ over the 
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reaction z/g + p — * n + e^ 



2. Calculations 

The model of the expansion of matter occurring in the fast winds is based on the previous 
fast expansion calculations of Meyer (2002) and similarly uses the Clemson Nucleosynthesis 
Code (Meyer 1995) with NACRE (Angulo et al. 1999) and NON-SMOKER (Rauscher & 
Thielemann 2000) rate compilations. In these calculations, material expands at constant 
entropy from high temperature and density. In the present parameterization, a fluid element 
in the wind moves out homologously so that its velocity w oc r, and r grows exponentially 
in time with timescale r. This means the density declines exponentially with time on a 
timescale r/3. This approximates the acceleration phase of the wind (e.g., Qian & Woosley 
1996). 

We computed several fast expansions. The final overproductions are presented in Ta- 
ble 1. For all calculations, the initial composition of the starting material was a mixture of 
protons and neutrons which gave the desired initial value of Ye. This material was given an 
initial temperature of Tg = T/10^ K = 10.0 and an initial density consistent with the chosen 
value of entropy. The calculations were halted after the temperature had dropped below 
Tg = 0.01, by which time all capture reactions had long since ceased and nuclei were only 
beta-decaying back to stability. The values of r, s, and initial Yg are used when referencing 
specific calculations. 

We first focus on a particularly interesting expansion with s/ks = 145, r = 0.003 
seconds, and Y^ = 0.510, which produces a significant quantity of the light p-process nuclei 
^^Mo and ^^Ru. Figure 1 shows the nuclear abundances as a function of atomic number for 
four times in the expansion. A significant abundance of remarkably high-mass nuclei builds 
up early in the expansion. This surprising result arises from the fact that the abundances 
of the light nuclei ^H, ^H, and ^He are low at high entropies. As a result, nuclear flows 
assembling ^He, which proceed through these light nuclei, become too slow to maintain 
an equilibrium between the free nucleons and ^He. As the temperature falls, equilibrium 
between the free nucleons and ^He would increasingly drive the free nucleons to assemble 
into ^He; however, the slowness of the requisite reaction flows prevent this from occurring 
at the necessary rate. Though this equilibrium fails early (by Tg ~ 9), the system is still 
primarily composed of ^He. For this particular calculation, at Tg si 6 the system is roughly 
90% ^He, 6% protons and 4% neutrons by mass with only a slight dusting of heavy nuclei. 
By Tg ^ 5 the system is roughly 95% ^He, 3.5% protons, and 1.5% neutrons by mass. From 
Tg = 6 to Tg = 5, the number of heavy nuclei have increased from ~ 10^^ to 10~^ per 
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nucleon; thus, the heavy nuclei are still negligible in abundance compared to the protons, 
neutrons, and *^He. The critical point is that the disequilibrium between free nucleons and 
^He allows a large abundance of free nucleons to exist at temperatures that would normally 
see the free nucleons locked into ^He. The few heavy nuclei that have been forming at these 
times thus coexist with a large overabundance of free nucleons (Meyer 2002). Through a 
rapid sequence of neutron and proton captures these few heavy nuclei are driven to quite 
high nuclear mass, as shown in the Tg = 4.39 panel in Figure 1. The peak at Z ^ 58 arises 
from the fact that the nuclear flows dam up at the closed neutron shell at A^ = 82. 

As the temperature continues to fall, new heavy nuclei assemble from the abundant 
^He nuclei. The growing abundance of heavy nuclei is then increasingly able to catalyze the 
synthesis of ^He from free neutrons and protons through reaction cycles such as 
^^Ni(n,7)^'^Ni(n,7)^^Ni(p,7)^^Cu(p,Q;)^^Ni. Once the neutrons disappear (which helped hold 
the nuclear abundances at the higher mass, nonequilibrium distribution), the heavy nu- 
clei begin to disintegrate towards iron-group nuclei, the favored isotopes in the equilibrium 
appropriate at those conditions. The Tg = 3.93 panel of Figure 1 shows the abundance 
distribution shortly after this disintegration flow begins. The peak at Z = 50 arises from 
the fact that at this point the nuclear flow is damming up at the Z = 50 closed shell. By 
Tg = 3.80, the disintegration flow has broken through Z = 50 and is now dammed up at 
A^ = 50, and predominantly at the isotope ^^Mo. As the temperature continues to decline, 
proton captures deplete ^^Mo and create a large abundance of ^'^Ru and ^^Pd. These isotopes 
subsequently decay to ^^Mo and ^^Ru after reaction freezeout. 

It is essential to note the large flnal abundances of protons and ^He. These species, 
and even a residual quantity of neutrons, are abundant throughout the disintegration epoch. 
This high abundance of neutrons, protons, and ^He nuclei hinders the disintegrations, despite 
the high temperatures at the time the disintegrations are occurring (Tg ^4). As mentioned 
above, the high abundance of protons also leads to proton-capture reactions along the A^ = 50 
closed shell that modify the abundances at late times. 

Table 1 presents top ten overproduction values from six calculations of nucleosynthesis 
in fast, high-entropy, proton-rich expansions. Calculation (B) is the model discussed above 
(s = 145 kB, Ye = 0.510, and r = 0.003 seconds). As discussed above, the most overproduced 
isotopes are light p-nuclei, which are made as A^ = 50 progenitors. At the top of the list 
are ^^Ru (made as ^^Pd) and ^^Mo (made as ^^Ru) at overproduction factors of ~ 10^. Also 
signiflcantly produced are ^^Mo (made as ^^Rh) and ^^Mo, which, although depleted late in 
the expansion by proton captures, is ultimately produced off of the N=50 shell as ^^Ru. 

Calculations (A) and (F) show results for expansions with lower entropy (s = 140 ks) or 
lower Ye {Ye = 0.505) than in our reference calculation, (B). In both cases, the disintegration 



-5- 



flow began earlier in the expansion than in calculation (B) and allowed the nuclei to return to 
the iron-group isotopes before freezeout. On the other hand, the larger entropy in calculation 
(C) delayed the onset of disintegration and thereby strongly overproduced higher-mass p- 
process isotopes. In calculation (D), still heavier proton-rich isotopes formed. Rare isotopes 
such as ^^°Ta are particularly strongly produced in this expansion. Finally, in calculation 
(E), the entropy was sufficiently high that the neutrons disappeared only at low temperature 
so the disintegration flow never occurred. This proton-rich expansion was able to produce 
heavy, neutron-rich isotopes whose production is normally associated with rapid neutron 
capture nucleosynthesis (cf. Meyer 2002). These results show the remarkable sensitivity of 
the nucleosynthesis to small changes in the expansion parameters. 

Because of this sensitivity and the non-equilibrium aspect of the nuclear flows during 
key epochs of the expansion, we expect that, for a given set of expansion parameters, the 
nucleosynthetic yields of such proton-rich expansions will be particularly sensitive to nuclear 
reaction rates on a large suite of isotopes. This is confirmed in Table 2, which show results 
for calculations identical to calculation (C) but with all charged-particle and electromagnetic 
rates on nuclei with Z > 27 increased or decreased by a factor of two. In the former case, 
the larger reaction cross sections enhance the assembly of nucleons into "^He and cause earlier 
and more efficient disintegration of the heavy isotopes into iron-group nuclei. In the latter 
case, the assembly of nucleons and, hence, the disintegration flows are less efficient and a 
heavier distribution of nuclei results than in calculation (C). Reaction surveys will be needed 
to further clarify this issue. 



3. Implications 

As Table 1 shows, there is a considerable variety of final abundance yields arising in 
these expansions. We do not expect these fast, proton- rich expansions to be dominant 
contributors to the solar system's r-process isotopes: the overproductions are low and the 
final abundance pattern is different from the solar system's r-process pattern. On the other 
hand, these expansions can make interesting quantities of rare isotopes (such as ^^"^Ta) and 
p-process nuclei, perhaps most intriguingly, the light p-process isotopes. 

The precise mechanism for the production of the light p-process nuclei, ^^'^^Mo and 
^^'^^Ru, remains mysterious. While the heavier p-process isotopes are well accounted for by 
the "gamma-process" in either core-collapse supernovae (Woosley & Howard 1978) or perhaps 
in thermonuclear explosions (Howard et al. 1991), the same nucleosynthesis underproduces 
the light p-isotopes. This is because, unlike the heavy p-process nuclei, the light p-process 
isotopes are nearly as abundant as their r-process and s-process counterparts, which serve as 



-6- 



seeds for the gamma-process. This suggests some other process than the "gamma-process" 
may be responsible for their origin (ahhough questions remain about the need for an exotic 
site (Costa et al. 2000)). Suggested other sites or processes are thermonuclear supernovae 
(Howard et al. 1991), alpha-rich freezeouts in core-collapse supernovae (Fuller & Meyer 1995; 
Hoffman et al. 1996), or the rp-process (e.g., Schatz et al. 2001). Each of these processes or 
sites has its own difficulties in either producing the right isotopes or actually ejecting them 
into interstellar space, and the question of the site of origin of the light p-process nuclei 
remains. 

As evident from Table 1, the typical overproduction factors for calculations that produce 
light p-process isotopes are ~ 10^. Core-collapse supernovae are responsible for production 
of ^^O, and such supernovae typically overproduce this isotope by a factor of about 10 
(e.g., Woosley & Weaver 1995). If we assume light p-isotopes are 1) produced in rapid, 
high-entropy, proton-rich expansions in all core-collapse supernovae, 2) are overproduced 
at a level of ~ 10^, and 3) are diluted into 10 Mq of ejecta, then we find each supernova 
must eject ~ 10^'^ Mq of this high-entropy, proton-rich matter. This is comparable to the 
estimates of the total mass ejected in proto-neutron star winds (e.g., Thompson et al. 2001); 
therefore, such winds may contribute interesting amounts of light p-process nuclei to the 
solar system if they indeed eject mildly proton-rich matter with the right timescales and 
entropies. On the other hand, some other as yet uncharacterized astrophysical site may 
achieve the necessary conditions. A full mapping of parameter space is clearly needed to 
understand this nucleosynthesis process as well as an exploration of its possible astrophysical 
settings. 

Because winds from proto-neutron stars stars may evolve from slightly proton rich to 
neutron rich with time, p-process and r-process nucleosynthesis may occur in the same site 
with the former preceding the latter by a few tenths of a second, the timescale on which 
the neutrino and antineutrino spectra are changing. This may have interesting implications 
for the coupling or decoupling of r-process and p-process isotopes in cosmochemical samples 
(e.g., Yin et al. 2002). 

Finally, we note that comparison of observations of abundances of Sr, Y, and Zr in 
low-metallicity stars and models of the chemical evolution of the Galaxy hint at a "lighter 
element primary process" (LEPP) that may contribute to the synthesis of A ^ 90 nuclei 
(Travaglio et al. 2004). The process we describe is primary and could therefore be an 
interesting component of this inferred LEPP. 

The authors are grateful to D. D. Clayton, M. D. Leising, and J. R. King for discussions 
and to L.-S. The for assistance. This work was supported by NSF grant AST 98-19877, 
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Fig. 1. — Elemental abundance versus mass number A for four times during the r = 0.003 
s, s = 145 ks, Ye = 0.510 expansion. 
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Table 1. Ten most overproduced isotopes in calculations of fast, high-entropy, proton- rich 

expansions 



rank 



O 



O 



O 



O 



O 



O 



1 


59 Co 


5.31 


lOi 


94 Mo 


2.15 


106 


"2Sn 


7.91 


106 


l«'Ta 


7.71 


10^ 


201 Hg 


7.91 


105 


59 Co 


2.78 


102 


2 


45 Sc 


4.98 


lOi 


96Ru 


7.74 


105 


ll3ln 


3.88 


106 


181 Ta 


1.27 


106 


200 Hg 


5.34 


105 


60Ni 


3.95 


IQi 


3 


60Ni 


4.25 


lOi 


95Mo 


3.87 


105 


108 cd 


2.09 


106 


176 Lu 


7.70 


105 


204Hg 


4.28 


105 


45 Sc 


2.59 


IQi 


4 


49 Ti 


3.40 


lOi 


92Mo 


5.40 


104 


iiOCd 


1.03 


10" 


175 Lu 


7.44 


105 


199Hg 


4.25 


105 


49Ti 


1.94 


IQi 


5 


63Cu 


1.64 


lOi 


93Nb 


2.06 


104 


106 Cd 


8.52 


105 


174Yb 


5.31 


105 


203 -pj 


3.25 


105 


63 Cu 


1.84 


IQi 


6 


48 Ti 


1.31 


lOi 


84 Sr 


1.22 


104 


i26Xe 


4.86 


105 


178Hf 


5.23 


105 


202Hg 


3.18 


105 


48Ti 


1.46 


IQi 


7 


44 Ca 


1.09 


lOi 


97Mo 


1.17 


104 


"4sn 


4.33 


105 


171 Yb 


5.22 


105 


198pt 


1.94 


105 


44Ca 


1.17 


IQi 


8 


43Ca 


8.28 


10" 


98Ru 


6.47 


10^ 


"5sn 


3.27 


105 


182^ 


4.15 


105 


182w 


1.74 


105 


43 Ca 


1.09 


IQi 


9 


47Ti 


6.13 


10" 


91 Zr 


4.39 


103 


124Xc 


1.59 


105 


184^ 


3.61 


105 


184w 


1.61 


105 


47tj 


6.41 


lOO 


10 


50 Cr 


6.11 


10« 


90 Zr 


2.37 


103 


102pd 


1.26 


105 


176Hf 


3.54 


105 


179Hf 


1.55 


105 


Sly 


6.08 


lOO 



'^Calculation A: t = 0.003s, s=140, initial Ye=0.510 
''Calculation B: r = 0.003s, s=145, initial Ye=0.510 
^^Calculation C: r = 0.003s, s=150, initial Ye=0.510 
^Calculation D: r = 0.003s, s=160, initial Ye=0.510 
^Calculation E: r = 0.003s, s=170, initial Ye=0.510 
^Calculation F: r = 0.003s, s=145, initial Ye=0.505 
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Table 2. Ten most overproduced isotopes in calculations of with modified reaction rates. 



ank 


^Z 


O 


^Z 


O 


1 


«Sc 


4.19 


lOi 


126Xc 


1.66 


10^ 


2 


60Ni 


2.96 


lOi 


I24xe 


7.62 


106 


3 


49 Ti 


2.92 


10^ 


"2Sn 


4.29 


106 


4 


59 Co 


1.38 


lOi 


l32Ba 


3.53 


106 


5 


63 Cu 


1.27 


lOi 


l^ORa 


3.53 


106 


6 


48 Ti 


1.20 


lOi 


ll^In 


2.44 


106 


7 


44 Ca 


9.89 


10" 


138 Cc 


1.96 


106 


8 


«Ca 


7.42 


10" 


l28Xc 


1.15 


106 


9 


47Ti 


5.41 


10" 


"4Sn 


9.30 


10^ 


10 


50 cr 


4.79 


10" 


"5Sn 


9.07 


10^ 



^Calculation A: t = 0.003s, s=150, initial Ye=0.510; charged particle and EM rates for Z > 27 X 2 
''Calculation B: r = 0.003s, s=150, initial Ye=0.510; charged particle and EM rates for Z > 27 X 0.5 



